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ABSTRACT 

In order to soften the effects of rub, the smooth 
stators of turbine gas seals are sometimes replaced by 
a honeycomb surface. This deliberately roughened 
stator and smooth rotor combination retards the seal 
leakage and may lead to enhanced rotor stability. 
However, many factors determine the rotordynamic co- 
efficients and little is known as to the effectiveness 
of these "honeycomb seals" under various changes in 
the independent seal parameters. This analysis 
develops am analytical-computational method to solve 
for the rotordynamic coefficients of this type of 
compressible-flow seal. 

The governing equations for surface roughened 
tapered annular gas seals are based on a modified Hirs' 
turbulent bulk flew model. A perturbation analysis is 
employed to develop zeroth and' first-order perturba- 
tion equations. These equations are numerically 
integrated to solve for the leakage, pressure, density, 
and velocity for small motion of the shaft about the 
centered position. The resulting pressure distribution 
is then integrated to find the corresponding rotor- 
dynamic coefficients. Finally, an example case is used 
to demonstrate the effect of changing from a smooth to 
a rough stator while varying the seal length, taper, 
preswirl, and clearance ratio. 
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H(z,e,t) 
h = H/C 
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damping coefficients defined by 
Eq. (26) 

= Shaft diameter 
= Local seal clearance 
= Dimensionless clearance 
= Entrance loss coefficient 


K, k 

K, k 

L 

l = L/R 


= Direct and cross-coupled stiffness 
coefficients of Eq. (1) 

= Dimensionless direct and cross-coupled 
stiffness coefficients defined by 
Eq. (25) 

= Seal length 

= Dimensionless seal length 



Pa 

(Rw) \ 


= Mach number 

= Coefficients for Hirs' turbulent 
lubrication equations 


= Pressure coefficient 


NOMENCLATURE 
C (Z) 

c = (CMC X ) n 

C * C/R 

c v 

C, a 

Ci a 


= Centered position seal clearance 
= Nominal seal clearance 
= Dimensionless nominal seal clearance 
= Specific heat at constant volume 

= Direct and cross-coupled damping 
coefficients of Eq. (1) 

= Dimensionless direct and cross-coupled 
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P 

P = P/P a 
R 


= Pressure 

= Dimensionless pressure 
= Shaft radius 


2pU C . . 

r = — = Centered position, axial Reynolds 

ao K number 

R c = p(Rw)H/u = circumferential Reynolds Number 

R = p (Rw) C/u = Centered [position, nominal 

circumferential Reynolds number 

R = Perfect gas constant 

T = Temperature 

t = Time 
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U = Rw 


u 0 = U 0 / (Rw) 
u z = U z /(Rw) 
Z, R9 


z = Z/L, 0 

V = C p /C v 
e 

e/2C 

V 
P 

P = P/P a 

V T r 

T =■ TW 

fi 

52 = 52/w 

U 

Subscripts: 
a, e, x, b 


0, 1 


s, r 


= Velocity of rotor surface 

= Fluid velocity in the z and 0 
directions 

= Dimensionless tangential and axial 
velocities 

= Axial and circumferential seal 
coordinates illustrated in Fig. (1) 

= Dimensionless seal coordinates 
= Specific heat ratio 
= Dimensionless seal eccentricity ratio 
= Relative surface roughness 
= Viscosity 
= Density 

= D intensionless density 
= Shear stress illustrated in Fig. (2) 

= Dimensionless time 
= Shaft orbital velocity 
= Shaft whirl ratio 
= Shaft angular velocity 


= Reservoir, entrance, exit, and sump 
conditions, respectively 

= Zeroth and first-order perturbations 
in the dependent variables 

= Stator surface and rotor surface, 
respectively 


INTRODUCTION 

Figure 1 illustrates the basic geometry of the 
convergent tapered annular turbine gas seal. In this 
figure both the rotor and the stator elements of the 
seal are shown to have the same nominally smooth sur- 
faces. In practice, however, the smooth stator is 
sometimes replaced by a honey comb or other deliberate- 
ly roughened surface. The purpose of this roughened 
surface is to soften the effects of rub from the rotor 
and to retard leakage. But in addition, the smooth 
rotor and rough stator ccmbination may have significant 
influence on the seal's rotordynamic coefficients. In 
fact, von Pragenau [1] suggests just such a concept for 
an incatpressible-flow "damper seal" which he believes 
will enhance rotor stability. 

As related to rotordynamics, seal analysis has the 
objective of determining the reaction force acting on 
the rotor as a result of the shaft motion. For small 
motion about a centered position, the relation between 
the -reaction-force c opponents and the shaft motion can 
be written 



The off-diagonal coefficients (k, a, m) are referred to 
as the cross-coupled stiffness, clamping, and added-mass 
terms, respectively. These cross-coupled terms arise 
fron the fluid's circumferential velocity component. . 
This phenomenon is usually referred to as the effects ol 



swirl. The circumferential velocity canponent is in 
part a function of the stator and rotor surface 
roughness . A rough stator and smooth rotor will tend 
to reduce the circumferential velocity, leading to a 
reduction in the destabilizing cross-coupled cerms. 

Fleming 12,3] made a separate analysis for the 
direct stiffness K, and for the direct danping C of 
smooth tapered annular gas seals. Hcwever, he did not 
include the effects of swirl and thus could not obtain 
the cross-coupled terms. Childs [4] developed an 
analysis for both direct and cross-coupled terms of in- 
ccmpressible-flow by using Hirs' [5] turbulent bulk- 
flow model and a perturbation technique. Nelson [6] 
used a similar approach to develop a numerical solution 
for the direct and cross-coupled stiffness and damping 
of smooth cctipressible-ficw seals. The present analy- 
sis modifies the solution of reference [6] to include 
the effects of different stator and rotor surface con- 
ditions and then demonstrates the analysis on a 
specific seal example. 


GOVERNING EQUATIONS 

The control volume element shown in Fig. 1 has 
been enlarged and redrawn in Fig. 2. Note that the 
smooth stator surface has been replaced by a roughened 
honeycomb surface. The shear stresses t s and x r are 
the net wall shear stresses resulting frcm both the 
pressure induced flow and the drag induced flow. Hirs' 
turbulent bulk-flow model assumes that these stresses 
can be written as 


T = *5PUV 
m o 



( 2 ) 


where is the mean flew velocity relative to the 
surface upon which the shear stress is acting. Hirs' 
constants Hq and n^, are generally snpirically deter- 
mined frcm pressure flow experiments. For the control 
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Figure 2. Control Volume for a Seal with 
a Honeycomb Stator. 
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Assuming a perfect gas (Cyl^p/p (y-1) and using the 
dimensionless parameters defined in the Nomenclature, 
the above governing equations take the following 
dimensionless working form. 

Momentum: 

-P_ u r "I Du 

— &=ZT f s + f r + ~ (4a) 

pi 3z ch 3 J Dt 


i e f s + ( V 


volume in Fig. 2, rig and mg represent Hirs' constants 
relative to the stator surface and n r and ity. represent 
those relative to the rotor surface. Substituting the 
mean flow velocity relative to each surface into 
Eq. (2) and then taking the appropriate component of 
the shear stress in the Z and 0-direction, the axial 
and circumferential momentum equations are 
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Ihe bulk-flow continuity equation is 

n 3(p H ) 1 3(PH 0 H > 3(pU z H) 

0 ~ 3t R 30 + 3Z 


f = — 
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f = — 
r 2 


«-2) 2 

(U 8 -1) 2 + U 2 2 


m r 

(2RC) r 


Dt 3t + w 0 30 + i 3z 


PERTURBATICN ANALYSIS 


Ihe governing Eqs. (4a) through (4d) define the 
relationship between the dimensionless clearance, pres- 
sure, density, axial velocity, and circumferential 
velocity (h, p, p, u z , uq) as functions of the inde- 
pendent dimensionless spatial variables (e, z) and the 
dimensionless time t'. Expansion of these equations 
in the perturbation variables 


h = h + eh. p = p + ep, 
o 1 o 1 


P = P + EP, u_ = u„ + EM, 


"0 = u 0o + “ei 


And for adiabatic flew the energy equation is 


yields the zeroth and first-order equations as shown 
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in Appendix A. 

Zeroth-Ordar Solution 

The zeroth-order aquations describe the steady 
flew resulting from a centered position rotating shaft. 
Before these equations can be integrated, values for 
Hirs' constants nig, ng, m,., and nj. must be established. 
Lacking experimental data for these constants, values 
can be approximated by the use of Colebrook's 
formula [7] . 



For a giyen relative roughness a least-squares fit is 
used to determine n 0 and 1 % over a range of % (say 
5000<R<1000000) . 

Integration begins by guessing an entrance zeroth- 
order Mach number M o (0) . Defining an entrance loss 
coefficient k in a manner similar to Zuk [8] , the 
following equations give the initial zeroth-order 
pressure and density. 


r (y-1) U+k)M 2 (of 

li + 2 ° J 

Y/(Y-1) 

1 + 

r (Y -i)M o 2 (o) 


L 2 


r (Y-l)Cl+k)M 2 (of 

Y/(Y-1) 

L 14 "— 

2 



In the first application, k is assumed to be 0.1. Ex- 
panding the Mach number as defined in the nomenclature 
in terms of the perturbation variables gives the 
following zeroth-order entrance equation 


along the seal length. Ihe guess for the entrance Mach 
number is continually adjusted until: (a) the Mach 
number at the exit just reaches one and the exit pres- 
sure remains greater than the sump pressure for choked 
flew, or (b) until the exit pressure just matches the 
sump pressure and the Mach number remains less than 
one for unchoked flow. All intermediate values of the 
pressure, density, and velocities and their deriva- 
tives are then stored for later use in solving the 
first-order perturbation equations. Also, the leakage 
is determined from these zeroth-order values. 

First-Order Solution 

Ihe first-order Eqs. (A. 2a) through (A. 2d) define 
PI (2/6, t), Pl (2,8/t) , Uj.i(z,e,T), and u^i (z,e,T> 
resulting from the seal clearance function (z,8,t) . 

If the shaft center moves in an elliptical orbit, then 
the rotation displacement vector to the shaft center 
has coordinates 

X = Cx Q cos.Qt , Y = Cy Q sinSlt (11) 

and the clearance function is 

h. = -x cosflxcosS - y sinPxsinS (12) 

lo o 

The assumed harmonic response is 

p. = (p c cx3sQT+p s sinflT)cos8+(p°cosnT+p s sinfiT) sin8 

1 x x y y 

p. = £cosnx+p s sin£2x ) cos0+ (p^cos^x+p^sinfix) sin9 
x x x y y 

u zl = (u^cosSx-m^sinSx) cos6+ (u^cosfix+UySinQx) sin6 

u„. = (v c cos!2x+v s sinfix)cos0+ (v c cosfix+v s sinI2x) sinO 
si x x y y 

(13) 

Substitution of Eqs. (12) and (13) into the first-order 
Eqs. (A. 2a) through (A. 2d) yields sixteen first-order 
ordinary differential equations which can be written 
in the form 


dx 

[A(z)J ^ + [B(z)]X = x Q C(z) + y^tz) (14) 


m o 2 (°) 


U^WgptO) 

YP 0 ( °)P C 


(9) 


Frcm this equation, the initial zeroth-order axial 
velocity u^lO) can be found. Having new p o (0), p o (0), 
and u zo (0) , the centered position axial Reynolds 
number R ao is determined and used to approximate a new 
loss coefficient from the data of Deissler [9] as 

k = >j5.3/log 1Q R ao - 1.0 (10) 

For Rao > 200,000, k is set equal to zero. This new loss 
coefficient is then used to determine new entrance 
conditions and the procedure repeated until a consis- 
tent result for k is found. Finally, the initial 
zeroth-order circumferential velocity, uq( 0), is a 
given independent variable which indicates the amount 
of prerotation given to the entering fluid. 


where 


„ ,~c ~s ~c ~s ~c ~s ~c ~s c s c s c s c s. 
X- (P X .P X /P y >Py/ P X # P x / Py/ Py'VVVVVVVV ‘ 


The coefficients of [A] , [B] , C, and D are given in 
Appendix B. These coefficients are completely deter- 
mined frcm the values obtained in the zeroth-order 
solution. 


Ihe necessary sixteen boundary conditions for 
Eq. (14) are new written by examining the perturbation 
conditions that must exist at the entrance and exit 
for choked or unchoked flow. For ease in writing 
these conditions, the following definitions are made: 
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(15) 


Having new the zeroth-order initial conditions , 
Eqs. (A. la) through (A. Id) are numerically integrated 


344 



(cLl) For choked flew, the first-order perturbation in 
in the exit Mach number is set equal to zero. 
Using the definition of the Mach number, the 
first-order perturbation is 


Pi u . 
M^ = M / — + — 


2P„ 


zo 


This yields 



(16) 


(a. 2) 


(b) 


V 1 * 

2p o (1) 


Hzld) 

u zo (1) 


h ll) 

2p 0 (i) 


(17) 


For unchoked flow, the first-order perturbation 
in the exit pressure is zero giving 

p^l) = 0 (18) 


At the entrance, the circumferential velocity 
perturbation is zero, i.e. 


u 01 (O) = 0 


(19) 


NUMERICAL EXAMPLE 

For the compressible flow seal with different 
stator-rotor surface roughness treatments, there are 
fourteen independent geometric and fluid dynamic seal 
variables. Assuming a perfect gas, these variables 
can be reduced to ten dimensionless parameters. One 
possible set is P 0 =p a / (p_w2R2>, R^ a = pgCuR/u, 

P a /Pb ' y f i>/D , C^/Cg, C/R, e s / 2C, e^/S? and 
uq(0)=Uq(0)/(wR) . Due to this large number of inde- 
pendent parameters, it is unreasonable to attempt to 
describe the leakage and dynamic coefficient's depen- 
dence in the form of a complete set of design charts. 
Thus, a specific seed geometry and flew condition was 
chosen, and only the length, taper, fluid prergtation, 
and clearance ratio (L/D , Cg/C*, ug o (0), and C/R) 
were independently varied. The particular seal 
selected is equivalent to the turbine interstage seal 
of the High Pressure Oxidizer Turbopump (HPOTP) of the 
Space Shuttle Main Engine operating .at Rated Power 
Level. The rotor is smooth and the stator is a honey- 
comb surface, resulting in the following seal 
parameters as supplied by Jackson [10] : 


(c) Expansion of the pressure loss Eq. (7) in terms 
of the perturbation pressure and the perturba- 
tion Mach number from Eq. (16) yields the 
following first-order pressure loss equation 
which must be satisfied at the entrance: 

- ... E P o (0) (h {0) , ^zl* 0 * 

p (0) + I + 

2p o (0)-r \P o (0) u 2Q (0) 



p = = 34.05 MPa (4938 psia) 
p^ = 26.41 MPa (3830 psia) 
T = 773°K (1391°R) 

a 

R = 7.282 an (2.867 in) 
C g = 0.381 mm (0.015 in) 


where 


Y (k+l)p Q (0 )M q (0) 


1+(Y-1) (k+1) 


m o 2 (°) 


( 21 ) 


(d) A similar expansion for the density change at 
the entrance defined by Eq. (8). gives 


^(0) + 


rsp Q (0) 


2u ?1 (0 ) Ei (°> \ 


2p o (0)+rs \ Uzo (0) p Q (0) 


where 


s = 1 +(y-1) 


[M q (0) (k+1) -2] 


= 0 
( 22 ) 

(23) 


2y (k+1) 

Solution of the differential Eq. (14) in terms of the 
above sixteen boundary conditions can be found by 
numerical integration techniques. The solution will 
take the form 

X = x^f (z) + y^tz) (24) 

Dynamic Coefficient Definitions 


As shown in reference [6] if the added mass terms 
are neglected, the dynamic seal coefficients can be 
obtained by numerically integrating the appropriate 
solution component of Eq. (24) , i.e. , 


1 1 


KC f 

P RL “ * J 

f^(z)dz 

£ _ kC _ 

K = = “TT 

P RL 

f f 3 (z;dz 

a 0 


a 

0 (25) 

1 



1 

cc a _ f 
p RL * J 

g 3 (z)dz 

•t=. 

€ 

11 

J" g 3 (z)dz 

0 


-a 

0 (26) 


C x = 0.254 mm (0. 010 in) 

L = 2.527 cm (0.995 in) 

Y = 1.4 

R g = 2480 m-N/kg°K (461 ft-lb/lbm°R) 

U = 2.05 x 10~ 5 Pa-s (1.38 x 10 -5 lb /ft-s) 

m 

U 9o (0) = 0.25 

(li = 28352 rpm 

e g /2C = 1.54 x 10" 2 =>m g - -0,0251, n g =0. 01534 

e /2C = 3.08 x 10 -4 =*m = -0.1691, n =0.03976 
r r r 

Results obtained for this seal are shown in 
Table 1. The first rew of data represents the seal 
with the given smooth rotor and rough stator. 'Ihe 
second row represents the results if both the stator 
and rotor are smooth. As expected, the rough stator 
decreases the leakage and cross-coupled stiffness. 
However, it also has the effect of significantly- 
reducing the direct stiffness and slightly decreasing 
the direct damping. 



P 

Kxl 0 -7 

feao -7 

0 

Q 


(kg/s) 

(N/m) 

(N/m) 

(N*s/m) 

(N-s/m) 

rough 

1.23 

3.07 

0.147 

1770 

13.0 

smooth 

1.35 

3.40 

0.160 

1300 

-0.9 


Table 1. Leakage and Rotordynamic 
Coefficients for the HPOTP Turbine Interstage Seal 


Results obtained by varying the seal geometry are 
plotted in Figs. 3 through 6. In these graphs, the 
broken lines represent the given rough stator and the 
continuous lines represent a smooth stator. The 


345 



Leakage and Rotordynamic Coefficients Leakage and Rotordynamic Coefficients 
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Fluid Prerotation, u, (0) 

00 


Figure 3. Leakage and rotordynamic coefficients 
vs. fluid prerotation (unchoked flow) 



Figure 5. Leakage and Rotordynamic Coefficients 
vs. Convergent Taper (unchoked flow) 
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Figure 6. Leakage and rotordynamic coefficients 
vs. length ratio (unchoked flow } 


Figure 4. Leakage and Rotordynamic Coefficients 
vs. Clearance Ratio, (unchoked flow) 
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vertical broken line represents the actual value of 
the independent variable for the HPOTP seal. Cross- 
coupled damping is not shewn since it was found bo be 
relatively insignificant. 

Figure 3 shows the effect of fluid prerotation 
(preswirl) . Clearly, prerotation has no effect on 
P, K, or C. However, there is a direct linear 
relationship for k. 

Figure 4 shows the effect of changing the naninal 
clearance C (convengent taper was held constant) . 
Within the range shown, decreasing the clearance 
results in an exponential increase in all coefficients 
and a linear decrease in leakage. 

Figure 5 shows the effect of convergent taper, 
Ce/Cx. For these results, the exit clearance was 
held constant and the entrance clearance increased. 

It should be noted that this also has the effect of 
increasing the naninal clearance. Thus 1. 0<C e /C > £_2. 5 
results in 0.0035<C/R<0.0061. As might be expected 
from Fig. 4, increasing the taper in this manner 
increases p and decreases k and C. However, K shews 
roughly a 50% increase when the seal is changed from 
straight to having a convergent taper ratio of 
Ce/Cx=2. 

Finally, Fig. 6 shows tire effect of seal length. 
Generally, as L/D increases, the coefficients increase 
and the leakage decreases. However, for a very long 
seal (i.e., L/D =0.8), k does reach a maximum and 
thereafter decreases. 

CONCLUDING REMARKS 


5. Hirs, G. G. , "Fundamentals of a Bulk-Flow Theory 
for Turbulent Lubrication Films," Ph.D. Disserta- 
tion, Delft Technical University, The Netherlands, 
July 1970. 

6. Nelson, C. C. , "Rotordynamic Coefficients for 
Compressible Flew in Tapered Annular Seals," 
Mechanical Engineering, Texas A&M University, 

1983. 

7. White, F. M. , Fluid Mechanics , McGraw-Hill, 1979. 

8. ' Zuk, J. , Ludwig, L. P., and Johnson, R. L. , "Quasi- 

One-Dimensional Compressible Flew Across Face 
Seals and Narrow Slots," NASA Technical Note 
D-6668, May 1972. 

9. Deissler, R. G. , "Analysis of Turbulent Heat 
Transfer and Flew in the Entrance Regions of smooth 
Passages," NACA IN 3016, 1953. 

10. Jackson, E. D. , Manager, Rotating Machinery Analy- 
sis, Rockwell International, Rocketdyne Division, 
personal correspondence, 1983. 

APPENDIX A: PERTURBATION EQUATIONS 


Zeroth-Order Eiquations 
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An analysis has been presented which calculates 
the leakage and rotordynamic coefficients for tapered 
annular gas seals m which the rotor and stator have 
different surface roughness treatments. To demon- 
strate this analysis, the effect of changes in seal 
length, taper, clearance and fluid prerotation was 
shown for the HPOTP turbine interstage seal. General- 
ly, changes in the abovementioned seal parameters 
resulted in major changes in the leakage and 
rotordynamic coefficients. 

In terms of the honey comb stator enhancing 
rotor stability, the results appear mixed. There is a 
favorable 9% reduction in cross-coupled stiffness and 
leakage. But at the same time, direct damping de- 
creases almost 2% and direct stiffness decreases 10%. 
Thus, general statements concerning the problems of 
instability and critical speeds can only be addressed 
by considering the entire rotordynamic system - clearly 
a problem outside the scope of this analysis. It 
should also be kept in mind that the selection of Hirs' 
constants may not accurately reflect the actual shear 
stresses developed over the honey ccmb surface. 
Experimental tests need yet to be performed to 
determine the correct values for these constants. 


3h 


3u 


3 Pr 


0 p o u zo 3z + p o h o 3z + u zo h o 3z 


(A.lc) 


3u 


0 = u 


6o 


. u . v fe, _ 

90 3z 20 3z 5 (Y-l) \*z 


* (u eo- 1)£ ro 

■Vzo 


3z 

(A. Id) 


where 


m +1 
s 


f__ = -sr 


f = — 
E ro 2 


2 2 

u + u 
o zo 


“s»» 


'V 


m r +l 


First-Order Equations 


(2R ) 
co 


references 

1. von Pragenau, G. L. , "Damping Seals for Turbo- 
machinery," NASA Technical Paper-1987, George C. 
Marshall Space Flight Center, Alabama 35812. 

2. Fleming, D. P. , "Stiffness of Straight and 
Tapered Annular Gas Path Seals," J. Lubr. 

Technol. , Vol. 101, No, 3, PP. 349-355, July 1979. 

3. Fleming, D. P. , "Damping in Ring Seals for Com- 
pressible Fluids," NASA CP-2133, 1980. 

4. Childs, D. W. , "Finite Length Solutions for the 
Rotordynamic Coefficients of Convergent-Turbulent 
Annular Seals," Mechanical Engineering, Texas A&M 
University, 1982. 


- p c 3 Pl 'Vi 9 P 0 


p Q 2 3z 


p*l 3z 


ch L 


u (f ,+f ,) 
zo si rl 


u ,-u t— 
zl zo h 


3u 


zl 


o 

dU' 


< f sl +f rl> 


. u 3u , u . 3u 
_ _ zl zo zl zl zo 

3x + u 6o 36 l 3z 2. 3z 


(A. 2a) 


347 



- p c 3 Pl 1 


P o 30 Ba o 


~ 'so) 
o 

+ <u 0o _1) ( f n-r f ro) +u 01 (f so +f ro> 


, 3u ei + u 3u ei u zo 3u 6i , u zi 3u eo 
+ ~ 8o “ + ~~ + “F~ 


(A. 2b) 


a = (m s +l)/(u 0o 2 + u zQ 2 ) 


b= (V 1)/t(u 0o- 1>2+u ao 2] 

APPENDIX B: MATRIX COEFFICIENTS 

All coefficients are zero except those defined 
below. Also, Y=YP C /[P 0 (Y“1)] and 
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(A. 2c) a i +12 ,i = Y ’ u zo/X. 

i=l,2,3,4 

3u . 
zl 

- Y 'Po u zo 

a i+8,i 

i=5,6,7,8 


3u_ 


3u, 


+ 2u 


zl 


+ u 


0o 


3z 


01 


3z 


3z 

3u, 


+ u 


0o 


01 


+ u 0o + ~ ( u ^ 
60 30 i V 20 3z 


L 3t 
3u 


0o 


J zl 


3z 


YP c 1 3 Pi , u !£i + ^0o!!i 
3 (Y-l) 1 Y 3 t 80 30 X. 3z 


, 3p 

+ 2 i u -u 

i 3z l 21 20 p 


. 3t 


tU8o !i + ^fi + i!!= 

8 30 l 3z l 3z ' 2l 



jcl _ ^11 
f ro h o 
(A. 2d) 


where 


= f 

L 1 

f h + h\ 

SO 

si 

V p o h o / 


f - = f 
rl ro 



r p i 

M 

m 

— + — 

r 

^ p 

h / 

L_ 

\ 0 

0 / 


+ a(u 0o u 91 + u zo u zl> 


+ b((u e o-i) u 6 i + u zo u z i> 


a = _S£ 
i+r,i i 


. _ U zo U 6o 

a i,i ~ l 


Matrix [B] 


b. 


1 d(U zo h o> 


i,i+4 «h Q dz 


i=9,10,ll,12 
i=13, 14,15,16 

i=l,2,3,4 


b i,i+5 _ -b i+l,i+4 = i-D3,5...11 

b i,i+6 = -b i+2,i+4 = u 0o 1=1 ' 2 ' 5 ' 6 ' 9 ' 10 
1 d(ph ) 


o o 


i=l,2,3,4 


b i,i+8 ~ Sh dz 
o 

b i,i+14 = -b i+2,i+12 = p o 1=1,2 

b i,i = ^T< m s f so^r f ro>- 37 ~ i = 5,6,7,S 

ch 0 p 0 Po 4 dz 




f so ll+ 4 a]+f ro (1 *4 b> 


j} + I 


du 

zo 

dz 


i=5,6,7,8 


b i,i+ 8 = ^ 2[,, 0o f so a+(u 0o- 1)f so b] i=5 ' 6 ' 7,8 

C) 


b i+6 , i - b i+8,i-2 ~ P c^ p o 


i=3,4 


b. 


-lu +{u -l)m f ] i=5,6,7,8 


1+4,1 dh 5 o s so o - r ro- 
o p o 

u 1 du. 

b i,i = =f C "eo f so a+(u 6o- 1)f ro b ) + I "dz 2 

o i=9,10,ll,12 
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”i,i« -^-I f a3 ‘ l «4»'«ro'WV- 1|2b ll 

i=9, 10, 11,12 


b i+12,i 


-YU do, 
zo 

p Q l dz 


h - H Y'fl 

i+ll,i ~ “°i+12,i-l y 

b i+10,i = -b i+13,i-l = Y ‘ u 


9o 


i=l,2,3,4 

i=2,4 

i=3,4 


b i+8,i 


Y ' U 2 o/ 2 P 0 dp O d Pnl 




P Q dz dz 


b i+7,i -b i+8,i-l = _Y * 


P o « 


p o 


r o 

i=5, 6,7,8 
i=6,8 


b = -b = -y' P ° Ue ° 

i+6,i b i+9,i-l Y ~ 

p o 


i=7,8 


b 



2ii du u Q du n 

zo zo . 6o 0o 

D i+4,i 

l \dz 

p Q dz ) 

*. dz l dz 


+ qf ro u zo 


1=9,10,11,12 

b i+3,i 

= " b i+4,i-l 

= u zo“ 

i=10,12 

b i+2,i 

= _b i+5,i-l 

= U ZO U 0O 

i=ll,12 

b i,i+l 

b i+l,i - 

"to 5 

i=13,15 

b i,i+2 

~ ~ b i+2,i = 


i=13 , 14 


Vectors C and D 


C 1 =d 4 


1 d(p o U zo> 


«h Q dz 


c 2 = -d 3 = -p o S/h o 
c 3 " ^2 " "^zo^o 


c 5 = ^ = ~i 
OT o 


c 9 " d 12 tu 0o (nl s“ 1)f so +(u eo _1) <V 1)f ro ] 


q(m r -l) 

c 13 = d 13 = ~h 
o 

P cPo S 


c 14 “ ^15 


P h 
o o 
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